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ABSTRACT 

Cosmological simulations of galaxy formation appear to show a 'two-phase' character with a rapid 
early phase at z>2 during which 'in-situ' stars are formed within the galaxy from infalling cold gas 
followed by an extended phase since z<3 during which 'ex-situ' stars are primarily accreted. In the 
latter phase massive systems grow considerably in mass and radius by accretion of smaller satellite 
stellar systems formed at quite early times (z > 3) outside of the virial radius of the forming central 
galaxy. These tentative conclusions are obtained from high resolution re-simulations of 39 individual 
galaxies in a full cosmological context with present-day virial halo masses ranging from 7 x 10 11 MqIi -1 
< M vir < 2.7 x IO^Mq/i -1 (h=0.72) and central galaxy masses between 4.5 x 10 10 Mq/i 1 < M* < 
3.6 x 10 11 Mi^hT 1 . The simulations include the effects of a uniform UV background, radiative cooling, 
star formation and energetic feedback from SNII. The importance of stellar accretion increases with 
galaxy mass and towards lower redshift. In our simulations lower mass galaxies (M* < 9 x 10 10 Af Q /i~ 1 ) 
accrete about 60 per cent of their present-day stellar mass. High mass galaxy (M* > 1.7 x 10 M^h ) 
assembly is dominated by accretion and merging with about 80 per cent of the stars added by the 
present-day. In general the simulated galaxies approximately double their mass since z=l. For massive 
systems this mass growth is not accompanied by significant star formation. The majority of the in-situ 
created stars is formed at z > 2, primarily out of cold gas flows. We recover the observational result 
of 'archaeological downsizing', where the most massive galaxies harbor the oldest stars. We find that 
this is not in contradiction with hierarchical structure formation. Most stars in the massive galaxies 
are formed early on in smaller structures, the galaxies themselves are assembled late. 
Subject headings: cosmology: theory - dark matter - galaxies: evolution - galaxies: formation - 
methods: numerical 
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1. INTRODUCTION 

Our understanding of galaxy formation has made great 
advances in the last two decades driven - primarily - 
by technological progress. Both ground and sky based 
measurements have allowed direct observation of various 
phases of galaxy formation and evolution over cosmic 
time with some detail ed information now available at 
redshi fts z > 2 (e . g. ISteidel et al.l Il999t iPettini et al l 
200"lt iGenzel et alj l2006t IForster Schreiber et al.l 120061: 



Truiillo et al.ll2007t iKrieh : et alj|2008t Ivan Dokkum et al 
2008 1 ; iMarchesini et al.l 120091: IForster Schreiber et al 
2009) . Simultaneousl y with a quite definite cosmological 
mode l (ACDM, e.g. ISpergel et al.ll2007l IKomatsu et al.1 
120101) having gained wide acceptance, we can, with in- 
creasing accuracy, simulate the evolution of galaxies 
from realistic initial conditions, with numerical reso- 
lution (in mass, space, and time) and physical mod- 
eling approaching the necessary degree of refinement 
(e.g. ISommer-Larsen et al.l I2003t ISpringel fe Hernquist 
2003 1 : Springell 120051: INagamine et all 120051: I Naab et al 



2007t IGovernato et alj 120071; iPiontek fe Steinmetzl 120091: 



Scannapieco et al., 2009; Sa wala et al. H2010l:IAgertz et al.l 
2010: iSchave et al.ll2010h 

The overall results are reassuring, with simula- 
tions and observations agreeing - in gross outline 
- as to t he time evolution of star /galaxy forma- 
tion (e.g. INagamine et al.1 120061: ISchave et all 12009ft 
as well as the global attributes of the galaxies 



distribution 



an d spatial organi 



such as lumi nosity 

zation (e.g. | Cen fe Ostrikerl Il999t IKauffmann et al 
119991 : ISomerville fe Primackj 119991 : I Springel et al.l 1200^ 



IKeres et al.l I2009al: IGuo et al.l l2010aTl . Understanding 
the development of the internal structures of galax- 
ies has been far more difficult to achieve with respect 
to the galactic stellar m ass fractions (e.g. IKeres et al.l 
I2009at IGuo et al. 2010 b;) as well as kinematics and mor- 
phologies (e.g. lAbadi et ; alj|2003t I Governato et"al~l l201CH: 
IFeldmann et al.ll2010D . ~ 

The terms with which we might usefully de- 
scribe such develo pment are still c ontroversial (e.g. 
Meza et al.1 120031: INaab et al.l l2007t IGovernato et al.1 



2007HPiontek fc Steinmetzll2009D . In a hierarchically or- 



ganized universe it has been natural to focus on over- 
dense 'lumps' of dark matter gas or stars and to fol- 
low tiie_jner^e^history of these lumps. A recent paper 
by iHopkins et al.l (|2009[) shows how useful this picture 
can be. But this is not the only de s cription of g alaxy 
formation. For exa mple IKeres et al.l (|2005l l2009bH and 
IDekel et al.1 ((2009b) have focused on how convergent cold 
streams of gas lead to early star bursts and the formatio n 
of the cores of mass ive galaxies. INaab et ail ([20071 12009T) . 
Uoung et~a l. ( 2009) and others have used high resolution 
hydro si mulations to explo re this phase in greater detail 
(see also iMeza et al.l ()2005h for the accretion histories of 
stellar halos of disk galaxies) . 

One fundamental and useful distinction is to examine if 
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a given star in the final galaxy was made (from gas) close 
to the center of the final system or, alternatively, near 
the center of some other, distant system and accreted in 
stellar form to the final galaxy. This distinction is use- 
ful, e.g. for understanding the s i ze evolution of massive 
galaxies (iKhochfar fc Silkil2006atlNaab et al.l[2007L l2009t 
Bezanson et al . 2009; N ipoti et al.ll2009at Hopkins et al.l 
2010t iFeldmann et al. I I2010D . In the simulations pre- 



sented here wc find that most stellar particles in massive 
galaxies are formed at high redshift either far inside the 
virial radius (< 3kpc) near the forming galaxy center or, 
alternatively in small systems outside the virial radius of 
the galaxy at a given cosmic time. We characterize the 
first category of stars as made 'in-situ' and the second as 
accreted or formed 'ex-situ'. In-situ stars are made (by 
definition) near to the galactic center over an extended 
time period. They are made from dissipative gas and, for 
massive systems, pro bably have relatively high metalicity 
(|Zolotov et al.|[2010l) . The peak rate of star formation for 
this category may be relatively early and in fact is very 
early (z sa 4) for the most massive systems. 

On the contrary, the accreted stars are typically made 
at quite early times as well, outside the virial radius, 
but added to the parent galaxy late in its evolution. 
They are added typically at radii larger than the ef- 
fective radius, r > r c g, and are expected to be metal 
poor, since the y originated in lo wer mass, lower metal- 
licity systems dNaab etaDHQOl). The ex-situ stars ac- 
crete via an energetically conservative process and their 
final binding energy is transferred to other phases (gas, 
stars, and dark matter) rather than simply radiated away 
(| Johansson et al.ll2009h . 

This alternative way of envisioning galaxy formation 
has many corollaries and makes many observed facts eas- 
ier to understand. In massive systems we expect consid- 
erable growth in mass and radius at late times but little 
star formation, with the late forming stellar envelopes 
comprised of stars which are typically older than the 
stars in the bulk of the galaxy. Further we find system- 
atic trends with galaxy mass. As one considers systems 
of lower mass, the in-situ component becomes increas- 
ingly dominant and the period of in-situ star formation 
is stretched out from being a small fraction of the Hubble 
time to a large fraction thereof. 

The paper is organized as follows. In section [2] we 
describe our simulations in detail, as there will be follow- 
up papers using this set of simulations. In addition we 
here discuss the conversion efficiency of gas into stars 
for our simulated galaxies. In section [3] we examine the 
dependence of the ratio of in-situ formed to accreted stars 
on the galaxy stellar mass along with its implications. 
We go on to analyze the half-mass radii of the different 
stellar components of our simulated galaxies in section 01 
Finally, in section [5] we summarize our findings. 

2. SIMULATIONS 

2.1. The large-scale dark matter simulation 

To find candidate dark matter halos for later refine- 
ment we performed a dark matter only simulation of 
a cosmological volume with a comoving side length of 
72Mpch~ 1 including 512 3 dark matter particles with in- 
dividual masses of m p = 2 x 10 8 Af Q /i _1 . The box is 
large enough to provide a representative piece of the 



universe and the mass resolution fine enough to allow 
us to reliably find dark matter halos with ~ 10 3 par- 
ticles being more massive than ~ lO n M0/i _1 . The 
initial con ditions were creat ed using GRAFIC1 and 
LINGERS (IBertschingerl[l995h . assuming a AC DM cos- 
mology with nearly scale-invariant initial adiabatic fluc- 
tuations. The cosmological pa rameters are based on the 
3-year results from WMAP (jSpergel et al.l I2007D with 
o- 8 =0.77, Q m =0.26, ft A =0.74, h = 0.72 (= H o =l00h 
kms _1 Mpc _1 ) and the initial slope of the power spec- 
trum is ro s =0.95. The initial conditions were then 
evolved from a redshift o f z ~ 43 to z = using 
GADGET-2 (|Springelll200"5l ) with a fixed comoving grav- 
itational softening length of 2.52kpc h. -1 . The simulation 
data was stored in 95 snapshots separated by Aa = 0.01 
beginning at a cosmological expansion factor of a=0.06 
(z«43). 

At z=0 we identify halos with a friends-of-friends al- 
gorithm and deter mine their centers using the shrinking 
sphere technique (|Power et~aLll200l . We then use the 
radius where the mean density drops below 200 times 
the critical density of the universe (r V i r = ^200) to mea- 
sure the halo mass therein (m V i r = TO2oo). This results 
in a complete halo catalogue (nhaios = 41313) for halos 
more massive than 2 x 10 10 Mq/i _1 which have proper- 
ties t ypical for this kind of simulation (see iMoster et al.1 
120101 for a detailed analysis of this simulation) . In brief, 
we show the dark matter halo mass function at z=0 and 
z=2 in Fig. [T] alon g with the analytical prediction from 
ISheth et allBool where we find small variations at the 
high mass end due to the limited boxsize. The corre- 
sponding distribution of the dimensionless spin parame- 
ter 



A' 



J 



V2r 



V r r v 



(1) 



,j vtr v c 1 vir 



defined bv lBullock et al.l (|2001ft . is shown in Fig. H Here 
J is the total angular momentum within r V i r and V c is the 
halo circular velocity V c 2 = Gm V i r /r V i r . The distribution 
of an gular momenta is consistent with previou s simula- 
tions puUoSitanilOOj |Vit^tSir2][lQQ2) and can 
be fitted with a log-normal distribution 



P(A') 



0.01 



exp 



ln 2 (A'/A(,) 
2a 2 



(2) 



with best-fit values X' — 0.038 and a = 0.58. 

2.2. Refined Simulations 

For the higher resolution re-simulations of individual 
halos we trace all dark matter particles that are closer 
than 2 x r2oo to the center of the halo at z=0. Following 
the halo back in time we include all particles in the trac- 
ing process which are within 2 x r2oo of the halo center at 
any given snapshot. This ensures that halo encounters 
during the assembly of the halo of interest are always 
resolved. We found this to be an efficient mechanism 
to reduce contamination with massive boundary parti- 
cles. The traced particles define the region for which we 
have to generate higher resolution initial conditions. For 
the cuboid enclosing this region we compute the short 
wavelength modes of the per turbation spectrum using 
GRAFIC2 (|Bertschingenl200"l . Based on the new spec- 
trum we replace the low resolution dark matter particles 
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TABLE 1 

Central Galaxies 



ID 




1 ~ 


m c 


1 ~ 


' ' Hns / I' 1 **, 
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2.73 
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11.6 


7.67 
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11.1 


9.59 


11.6 


8.31 


200865 
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716832 


0162 
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26.2 


2.64 


0.129 


10.8 


8.49 


11.2 


2.58 


134454 


1 06554 


494315 


0163 


689 


309 


25.3 


1.73 


0.150 


10.5 


9.11 


10.8 


4.75 


139297 


119486 


536504 


0175 


699 
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26.5 


1.29 


0.270 


11.3 


9.74 


11.8 


9.56 


127745 


117170 


530274 
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511 
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22.7 


1.71 


0.146 


10.3 


8.39 


10.6 


3.81 


1 03075 


98844 


405894 
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19.3 


1.42 
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10.8 


8.77 


11.2 


8.31 


102722 


99548 
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14.4 
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10.9 


9.71 


11.5 


9.26 
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8.15 


1 00251 
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0.300 


10.1 
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11.3 


8.31 


49561 


50348 


200794 


0443 
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16.6 


1.91 


0.277 


10.3 


6.55 


11.7 


8.31 


50289 


52800 


210493 
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11.7 


0.93 


0.361 


10.8 


10.0 


11.2 


8.79 


48521 


40463 


181178 


0549 


216 
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8.38 


0.450 


0.262 


10.7 


8.71 


11.4 


9.41 


39034 


39605 




0616 
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9.38 


0.455 


0.367 


10.6 


9.88 


11.0 


9.72 


34520 


37188 


147962 


0664 
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197 


7.48 


0.558 


0.343 


10.7 


9.06 


11.6 


9.41 


34393 


30862 


1 38039 


0721 


147 
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9.63 


0.629 
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7.07 


11.0 


6.69 


22910 


34776 


1 1 6680 


0763 


150 
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9.85 


0.177 


0.197 


10.4 


9.19 


10.8 


6.37 


25283 


34151 


119180 


0858 
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181 


10.3 


0.790 


0.355 


8.92 


5.49 


10.8 


6.69 


21022 


33759 


1 1 0365 


0908 
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9.67 
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8.84 


6.55 


10.8 


7.50 


19927 
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0948 
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1 7"? 

1 t O 
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u.ouo 


in fi 


y.oo 


11 2 


y . ou 


22627 


9^1 A7 


0447^ 
y^^i ( o 


OQ^Q 
uyoy 
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6.05 


0.399 


0.371 


10.1 


9.46 


10.5 


9.41 


23591 


23027 


94670 


0977 


94.4 


159 


4.55 


0.464 


0.496 


9.10 


7.21 


11.0 


8.63 


16592 


20916 


75143 


1017 


106 


166 


6.39 


0.736 


0.584 


10.0 


8.92 


11.5 


9.87 


21049 


20634 


83999 


1061 


103 


164 


5.18 


0.439 


0.335 


9.98 


8.72 


10.6 


8.15 


19196 


20400 


81076 


1071 


106 


166 


7.79 


0.610 


0.317 


9.66 


7.06 


10.9 


8.15 


18696 


24045 


84818 


1091 


112 


169 


7.53 


0.416 


0.280 


9.24 


5.37 


10.7 


6.20 


18487 


26210 


89119 


1167 


93.0 


159 


7.37 


0.659 


0.331 


9.32 


5.88 


11.0 


6.85 


15966 


22371 


75141 


1192 


78.0 


150 


4.36 


0.157 


0.442 


10.4 


9.54 


11.0 


9.56 


13041 


15792 


60404 


1196 


95.4 


160 


7.73 


0.99 


0.490 


9.23 


6.96 


11.4 


7.67 


16839 


20987 


75883 


1646 


71.3 
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5.38 


0.509 


0.480 


8.90 


6.23 


11.4 


8.31 


11143 


16557 


56264 


1859 


70.0 


144 


4.52 


0.340 


0.429 


9.82 


7.86 


11.3 


9.56 


12355 


16458 


56488 



Note. — all masses in units of 10 'h~ 1 Mq, timescales in Gyr. 
Qvirial mass, Qvirial radius in kpc/h, Utellar mass inside rio ,Qga,s mass inside rio,Q"atio of in-situ to ex-situ created stars, Qmean stellar 
age , @mean stellar age of in-situ created stars, 0mean stellar age of ex-situ created stars, [Jookback time where 50 per cent of the final 
stellar mass is in place, Qnumber of gas particles inside r200 , Qnumber of star particles inside r200 1 Qtotal number of particles inside r200- 
The horizontal bars indicate the separation into small, intermediate and high mass galaxies used throughout this paper 




Fig. 1.— Dark matter mass Function (solid) of the (72/i _1 Mpc) 3 Fig. 2.— Spin parameter distribution for the dark matter box. 

box at z =0 (red) and z=2 (b lue). The dashed lines show the pre- The green line shows the log-normal-fit with best fit values \' = 

diction of lSheth et all pOOlT ). 0.038 and a = 0.58. 
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with dark matter as well as gas particles at higher resolu- 
tion (f2f,=0.044, ridm=0.216). We only consider coherent 
regions within the cuboid that actually contain traced 
particles. Other regions as well as a 'safety margin' of 
IMpch -1 around the high-resolution cuboid are popu- 
lated with particles from the original initial conditions. 
To approximate the long range tidal forces, particles from 
the original simulation being further away from the cen- 
ter are merged, with the particle masses increasing as the 
square of the distance from the region of interest. By this 
and the inclusion of periodic boundaries tidal forces from 
distant regions are accurately included in the computa- 
tions. 

We obtain amoeba shaped initial conditions (see Fig. 
[3]) for which, on average, approximately 30% of the high 
resolution dark matter particles end up inside the v irial 
radius at re dshift z = (see IPower et al.l 120031 and 
in particular Uenkinsl 120101 for alternative ways of cre- 
ating high resolution initial conditions). The particle 
number in the boundary region is kept low enough to 
perform the simulations in a reasonable time. For ex- 
ample the most massive halo #0040 which has a to- 
tal mass m2oo of 2.7 x 10 13 Af Q /i _1 at z — took 
~ 23000 CPU-hours to evolve (3.8 x 10 6 high- resolution 
particles in dark matter and gas each). In the re- 
simulations the particles in the high resolution regions 
have a gas and star mass of m^, gas — 4.2 x 10 6 M Q /i _1 
(we spawn one star per gas particle) and a dark mat- 
ter mass of rridm = 2.5 x 1O 7 M /i -1 which is 8 smaller 
than in the original simulation. The comoving grav- 
itational softening length for the gas and star parti- 
cles is 400pch _1 and 890pch~ 1 for the high resolution 
dark matte r particles, sca led with the square root of the 
mass ratio (|Dehnenll2001[) . Compare d to some other re- 
cent cosmological zoom simulations (iScannapiecq et all 
20091: iGovernato et al l 120091: iPiontek fc Steinmetd 120091: 
Feldmann et al.l 1201 Of ) the resolution level of our simu- 
lations at M} la i « 10 12 Mq is slightly lower. But while 
these simulations are limited to a few halos in a small 
mass range we performed a significantly larger number of 
re-simulations of halos spanning a mass range of almost 
two orders of magnitude. The present-day properties of 
our re-simulated galaxies can be found in Table [TJ Fi- 
nally, we also performed a number of re-simulations at 
higher resolution, i.e. with particle masses 8 times lower 
and half the softening length. While increasing resolu- 
tion can slightly change the individual accretion histories 
of the galaxies, the global trends found in this paper re- 
main the same. 

2.3. Simulation details 

The simulations presented here have been per- 
formed using t he parallel TreeSPH code GADGET-2 
(|Springell 12005) which calculates the gas dynamics us- 
ing the Lagrangian Smoothe d Particle Hydro dynam- 
ics technique (SPH, see e.g. iMonaghanl I1992D . The 
code ensures the conservat ion of energy and entropy 
(|Springel &: He rnquist 2002) and includes star formation 
and cooling for a primordial composition of hydrogen 
and helium, where the cooling rates are computed un- 
der the assumption that the gas is optica lly thin and 
in ionization equilibrium (|Katz et al.l 119961 ). Addition- 
ally, the simulations include a spatially uniform redshift- 




FlG. 3. — Central region of the initial conditions for halo #0408 at 
z=43. The innermost region consists of the high-resolution gas and 
dark matter particles (red and blue). The green particles are dark 
matter particles taken from the original dark-matter-only run. The 
outermost dark matter particles have increasing mass depending on 
the distance, with sufficient resolution to represent the long range 
tidal forces. 

dep endent UV background radiation field with a modi- 
fied [HaardFj^JVIadau 1996 s pectrum, where reionization 
takes place at z ~ 6 (|Dave et al.lll999l ) and the UV back- 
ground radiation field peaks at z sa 2 — 3. For a recent 
detailed investigation on the effects of varying the back- 
grou nd radiation field on the evolution of galaxies, see 
e.g. (|Hambrick et al.H2010D . 

For the star formation and feedback prescription 
we use the self-regulated su pernova feedback model of 
ISpringel fc Hernquistl 120031 This models treats the 
ISM as a two-phase medium (|McKee fc Ostrikerl [l977t 
[Johansson & Efstathiou 2006|) where clouds of cold gas 
are embedded in the hot gas phase at pressure equilib- 
rium. Stars are allowed to form out of the cold gas phase 
if the local density exceeds a threshold value (n > n t h = 
0.205cm -3 ) which is calculated self-consistently in a way 
that the equation of state is continuous at the onset of 
star formation. In this model the star formation rate 
is set by dp*/dt = (1 — /3)p c /i*, where /3 is the mass 
fraction of massive stars (> 8M Q ), p c is the density of 
cold gas and is the star formation time scale set by 
t* = ^(n/nth) -1 ^ 2 ■ The massive short-lived stars heat 
up the surrounding gas with an energy input of 10 51 ergs 
due to supernovae. In order to prevent spurious star 
formation at high redshift we require an over-density of 
A > 55.7 for star formation to set in. 

2.4. The baryonic mass budget 

Using similar parameters for zoom simulations has 
been shown to result in galaxies with reasonable present 
day p roperties (|Naab et al.ll2007l l2009t [Johansson et all 
2009, 2010 in prep.). However, the employed star for- 
mation prescription favors efficient star formation at 
early times resulting in preferentially spheroidal systems 
with old stellar populations, due to the strongly self- 
regulated feedback which does not produce the super- 
nova driven winds. Fig. [4] shows the conversion ef- 
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Fig. 4. — Fraction of baryons that is converted into stars at red- 
shift zero. The vertical dotted lines indicate the separation into 
the different mass bins. The solid black li ne shows the results of 
IGuo et alJ !|2010bl) . the dashed line those of lMoster eTaLl $20Wt ). 

ficiency of the simulated galaxies at the present day 
/* = * rn V i ri dark) where m* is the stellar mass 

within 10 % of the virial radius, /b = Slb/f^m = 0.20 
is the cosmic baryon fraction and vn V i r ^ark is the dark 
matter mass within the virial radius of the galaxy. There- 
fore fb * m V i rt dark is the amount of total baryonic matter 
available in each halo and /* the fraction thereof that 
is converted into stars in the central galaxy. This frac- 
tion declines in a roughly linear fashion with the loga- 
rithm of the halo mass from /» w 0.5 for the smallest 
halos (« 10 11 - 9 M Q ) to R3 0.15 for high mass halos 
(> 10 13 M Q ) , still over- predicting by a factor of 2 the 
estimation from recent models (see however ? who find 
higher efficiencies for high mass galaxies) that are tested 
by matching observed l uminosity functions to simu- 
lated halo mass functions dMoster et al.ll2010tlGuo et all 
I2010ri IConrov fc Wechslerl 120091: IBehroozi et al l 2010) 
or weak lensing observations (|Mandelbaum et al. I f2006f) . 
Note that a Salpeter initial mass function would increase 
the 'observed' conversion efficiency by approximately a 
factor of two (?). 

The physical processes probably responsible for this 
discrepancy are well studied and it has been argued 
that feedb ack from SNII is important for low mass sys- 
tems (e.g. iLarsonl Il97l IDekel fc Silkl Il98l IGuo et all 
l2010a|) and feedback from super massive black holes 
dominates for h i gh m ass systems (jCroton et al.l 120061 : 
iDi Matteo et all 1200811 Although this issue is rela- 
tively well understood and many idealized calculations 
have shown how these feedback processes can expel 
the baryons from galaxies, there have been only a 
few high resolution galaxy formation calculations, us- 
ing cosmological initial conditions, beginning to mas- 
ter the physics well enough to match either the winds 
seen in forming galaxies or the fin al metal distribu- 
tion between galaxies an d the IGM ()Scannapieco et all 
120081: ISawala et all 120101 ) . Some other calculations do 
successfully allow for winds and for the consequences 
these winds have on the galaxies and the surrounding 
ISM (iSpringel fc Hernquistl 120031 : lOppenheimer fc Dave 
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Fig. 5. — Stellar origin diagram for all stars within rio at z = 
in galaxy M0125. Every grey dot indicates the time when a stellar 
particle was born and the distance, in units of the virial radius of 
the main galaxy at that time, where it was born. The contours 
enclose 90% (purple), 80% (turquoise), 60% (blue), 40% (green), 
25% (orange) and 10% (red) of the stars, respectively. The blue and 
red vertical lines show the virial radius and 10% of the virial radius, 
respectively. There is a clear distinction between stars initially 
formed in the galaxy and those formed outside the galaxy and are 
accreted later on (78 percent of all stars). 
(Dots excluded for file size reasons) 




Fig. 6. — Same as Fig. [5] but for the low mass galaxy M1646. 
There is significant in-situ star formation at the center even at low 
redshift and significantly less accretion of stars. In this case only 
52 per cent of the stars are accreted. 
(Dots excluded for file size reasons) 



Wiersma et all 120101: McCarthy et all 120101 ) . "Our com 



putations do not generate significant winds at high red- 



shift (e.g. iSteidel et al.ll20fol ) and thus overestimate, by 
roughly a factor o f two, the condensed baryon fraction o f 
massive galaxies (|Guo et alll2010bl: IMoster et alll2010D . 
This becomes worse if we extend the sample to lower 
masses where the halo occupation models predict a sharp 
drop off the conversion efficiency /* . This is probably due 
to the fact that ejective supernovae wind feedback, which 
is not included in the present study, is most effective 
in this regime. We are currently working to implement 
physically valid feedback implementations to address this 
problem. 

3. THE TWO PHASES OF GALAXY FORMATION 

The stellar particles ending up in the simulated galax- 
ies at z = are of two different origins. Some fraction 
of the stars are made in-situ, within the galaxies, from 
accreted gas and some fraction of the stars are made 
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Fig. 7. — Same as Fig. [5] but for all galaxies in low mass halos 
in the mass range 7.0 X 10 11 - 1.3 X lO 12 h _1 M (panel a), for 
intermediate halo masses in the range 1.3 X 10 12 -4.5 X 10 12 /i- 1 M Q 
(panel b), and for all high mass halos with 4.5 X 10 12 — 2.7 X 
10 13 /i — 1 Mq (panel c). The contours show the same percentiles as 
in Figs. OandE] The stars form in two phases, either inside rio or 
outside r v i r as can be seen in panel d. Gala xies in low mass halos 
have ongoing in-situ star formation (see Fig. HOC at relatively high 
specific rates until the present day, whereas in the highest mass 
group most star formation is complete by z=2. 




1 2 3 4 5 6 7 
Redshift 



Fig. 8. — Stellar mass assembly histories (solid lines) for low 
mass (top), intermediate mass (middle) and high mass (bottom) 
galaxies. The assembly is separated into in-situ stars (dotted line) 
and ex-situ stars that are accreted onto the galaxy later on (dashed 
line). The assembly of higher mass galaxies is more dominated by 
in-situ formation at high redshift, however, the total fraction of 
accreted stars by z=0 is higher (ss 80%) for massive systems than 
for low mass systems (ft 60%). 

ex-situ outside the galaxies and are acc reted later on 
(|Naab et al.ll2007t Uohansson et all 12009). The relative 
amount of in-situ and ex-situ stars is found to vary sys- 
tematically with galaxy mass. Two typical stellar origin 
diagrams indicating this behavior are shown in Fig. [5] 
and Fig. [6l 

To construct these diagrams we follow every star that 
ends up within 10% of the present-day virial radius of a 
simulated galaxy back in time. We use 10% of the virial 
radius, rio, as a fiducial value for the extent of the stel- 
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Fig. 9. — Average ratio of in-situ created stars that formed inside 
the bins indicated by the vertical dotted lines out of gas that was 
accreted cold to the total mass of in-situ created stars (dashed 
line) . The solid line shows the ratio of the stars created in-situ out 
of cold gas to the total stellar mass growth, this includes in-situ 
star formation as well as accretion. The error bars corespond to 
the lcr-dispersion. 



lar component of a simulated final galaxy inside its dark 
matter halo. Then we mark the time when a star was 
born as well as its distance from the galaxy center in units 
of the virial radius (at this time) with a grey dot. The 
values are discrete in time representing the discrete snap- 
shots. The contours in these plots encompass the small- 
est number of bins that include 10, 25, 40, 60, 80 and 90 
per cent of the stars, respectively. In FigjS] we show the 
stellar origin diagram for a massive system with a halo 
mass of ~ 10 13 M Q /i _1 . At redshifts z > 2 there are two 
separate peaks of star fomation: one inside rio, which is 
in-situ star formation and another one outside the virial 
radius of the system at that time. This indicates that a 
significant fraction of the stars in the present-day galaxy 
was made outside the galaxy and has been accreted later 
on. For this system the in-situ star formation decreases 
towards lower redshifts. Although there is ongoing star 
formation until z=0 the contribution to the final galaxy 
is negligible, since the contoured regions include 90 per 
cent of all stars in the galaxy. For a lower mass system 
with a halo mass of 7.1 x 10 11 M Q /i _1 the same analysis is 
shown in Fig. [6j In this case the fraction of stars forming 
ex-situ is lower and the contoured regions extend up to 
the present day, i.e. in-situ star formation continues at 
a significant level towards lower redshift. 

In Fig. [7j we have stacked all simulated galaxies of our 
sample into three mass bins (indicated by the horizontal 
bars in table [lj with the same number of objects (13), 
every star particle is weighted according to the total num- 
ber of stars in its host galaxy, so that every galaxy has an 
equal weight. The low mass bin contains galaxies with 
halo masses in the range 7.0 x 10 11 - 1.3 x 10 12 M Q /i _1 
(panel a), intermediate mass galaxies have 1.3 x 10 12 — 
4.5 x 10 12 M Q /i _1 (panel b) and high mass galaxies have 
4.5 x 10 12 - 2.7 x IO^Mq/i" 1 (panel c). These plots 
again demonstrate in a more statistical sense that the 
stars ending up in the final galaxies form in two distinct 
phases, namely in-situ in the galaxy and ex-situ outside 
the virial radii of the galaxies (red vertical dashed lines). 
The spatial division line between these two phases of star 
formation is at about 10% of the virial radius indicated 
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Fig. 10. — Star formation histories for low mass (upper panel), 
intermediate mass (middle panel) and high mass galaxies (lower 
panel) for all stars that end up inside the galaxy at z = 0. The 
solid black line shows the formation of the in-situ created stars, 
the red dashed line the formation of the ex-situ stars and the green 
dotted line shows the accretion rate of the ex-situ stars onto the 
galaxy. The arrows on top indicate the time at which half the stars 
are formed/added. The arrow at the bottom indicates the time at 
which 50% of the final galaxy mass is assembled. 

by the vertical blue dashed lines in Fig. [7] In addition, 
there is a clear trend that low mass galaxies have rela- 
tively more in-situ star formation at low redshift z < 1 
than higher mass galaxies. For the most massive galaxies 
the contribution from late in-situ star formation is rela- 
tively small. Panel d shows a histogram for the formation 
radii for all stars in all simulations. For this analysis we 
use 45 logaritmically evenly spaced bins. We see two 
peaks, for the in-situ created stars at log(r/r v i r ) « —2.5 
and for the ex-situ created stars at log(r/r v j r ) sa 0.6, re- 
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Fig. 11. — From top to bottom: Fraction of in-situ stellar mass vs 
total stellar mass inside rio, halo concentration and stellar density 
inside i?i/2 at Z = 0. There is a clear trend that galaxies with less 
in-situ star formation are more massive, have less concentrated 
halos and lower density central regions. 

spectively. A third peak appears between rio and r v ; r 
that is due to infalling substructure that is still star- 
forming. 

In Fig. [8] we show the average mass accretion histo- 
ries for the stellar particles in the three mass bins sep- 
arated into in-situ and ex-situ/accreted stars depending 
on whether they have formed inside or outside 10% of 
the virial radius. The galaxy growth is dominated for 
all three mass bins by in-situ star formation until z ss 2, 



when the mass of accreted stars equals the mass of in-situ 
stars. By z = about 41 ± 9% (we give mean values and 
the ler-dispersion of the 13 galaxies) of the stars in the 
low mass sample (top panel) have formed in-situ, the rest 
were accreted. For the intermediate mass galaxies (mid- 
dle panel) the fraction of in-situ stars is lower than for 
the low mass sample of ss 33 ± 10%, and 67% of the stars 
were accreted. With 78 ± 7% the fraction of accreted 
stars is even higher for the massive galaxies. On average 
only 22% of the present-day stellar mass is formed in-situ 
which is the dominant mode until z ~ 2 but thereafter 
contributes very little to the st ellar mass growth . 

Following iKereset all (|2005l ) and iKeres etaLl (|2009aft 
we examined whether the gas out of which the in-situ 
stars are formed in our galaxies was ever heated above 
Thot > 2.5 x 10 5 K throughout the simulation. The re- 
sults can be seen in Fig. [9j The dashed line shows, that 
up to redshift 2, where in-situ star formation is still dom- 
inating over accretion, almost all of the in-situ stars are 
formed out of gas that was accreted cold. Only at later 
times (0 < z < 2), when stellar accretion is the primary 
source of stellar mass growth, in-situ stars are forming 
out of cooling hot halo gas. At lower redshift the con- 
tribution of in-situ star formation out of cold flows to 
the total stellar mass growth becomes almost negligible 
(dotted line in Fig. [S]). The interpretation of the results 
does not change when we instead of a fixed temperature 
cut use a temperature thre shold related to th e current 
halo virial temperature ( see IKeres et al. ( 2005)). This is 
in agree ment with the previo us results of numerical sim- 
ulations (IKeres et al.l l2009af ) and analytical predictions 
(|Dekel fc Birnboimll2006D that galaxy growth at high red- 
shift (z > 2) is dominated by cold accretion. 

Fig. [TU] illustrates the star formation and assembly 
histories for the galaxies in the three mass bins. The 
red dashed line shows the archaeological star formation 
history of the accreted stars computed from the mass 
weighted ages of the accreted stars at the present day. 
All curves show a steep increase towards the peak at 
z « 4 at values of ~ 25M & yr~ 1 , « 55M Q yr~ 1 , and 
lbOMQ-yr^ 1 for the low, intermediate and massive 
bin, respectively. This is followed by an approximately 
exponential decline towards z = 0. The red arrow on 
top indcates the time when half of the accreted stars are 
formed. In all cases, i.e. at all masses this is at z rj 3. 
The green dotted line shows when these stars are ac- 
creted onto the galaxies. As this happens in mergers, the 
curves show peaks. On average the rates increase towards 
z — 2 and then stay relatively flat with average rates of 
Rj 3.6A/o2/r-\ « S^M^yr' 1 , and w nM^yr' 1 . The 
green arrow on top indicates when half of the present- 
day mass in ex-situ stars is accreted onto the galax- 
ies. This happens around z=0.7-1.2 and therefore signif- 
icantly later than the formation of these stars at z=3-4. 
The black solid line shows the formation history of the in- 
situ stars in the galaxies. This is most closely related to 
the star formation rate that would actually be observed 
in these galaxies. All curves peak at z > 3.5 at rates 
between m 5 and ss 20M & yr~ 1 . Independent of galaxy 
mass all rates drop to ss 2 — 3M0j/r _1 at z = 1 and stay 
constant to the presen t day similarly to th e observations 
of massive galaxies by Uuneau et al.1 (|2005l ). This results 
in a specific star formation rate of 0.31 ±0.15, 0.18 ±0.15 
and 0.053 ± 0.071 x lO^V" 1 for the different mass 
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Fig. 12. — Mean age of the stars inside rio as function of galaxy 
mass. High mass galaxies consist of older stars than the low mass 
galaxies, recovering the phenomenon usually referred to as 'archae- 
ological downsizing' (tmean oc log Mj' 6 ). 



bins. According to the definition by iFranx et al.l 120081 
(SFR/m* < 0.3/tthub) the galaxies in the high mass 
bin would correspond to quiescent galaxies. The time 
when half of the in-situ stars are formed is indicated by 
the top black arrows. This changes systematically with 
galaxy mass from z=1.4 to z=1.9 and z=2.1, i.e. the in- 
situ component is oldest for the most massive systems. 
The black arrow at the bottom of the panels indicates 
the time when half of the final galaxy was assembled. 
For all galaxies this is around redshift z « 1. Therefore 
all galaxies double their mass since then. For low mass 
systems the low redshift growth is dominated by in-situ 
formation whereas for high mass syst ems it is dominate d 
by accretion of small stellar systems (|Tiret et al.ll2010l ). 

In summary, at high redshift the assembly of galaxies 
at all masses is dominated by in-situ star formation fed 
by cold flows. The larger the galaxy mass the smaller 
is the late contribution of in-situ star formation. At low 
redshift, z < 1, the growth of low mass galaxies continues 
by in-situ star formation and stellar accretion whereas, 
massive galaxies gro w predominantly by accretion of ex- 
situ stars (see e.g. iFeldmann et al.l 120101 : iNaab et all 
f2009h . 

In Fig. [TT] we show interesting correlations of galaxy 
and halo properties with the fraction of in-situ stars indi- 
cating that this quantity is an important tracer of galaxy 
assembly. Essentially, this ratio, mi ns /m*, is a dimen- 
sionless measure for the degree to which the galaxy was 
formed by a dissipational v ersus a dissipationless process 
(jLackner fc Ostrikerll2010T ). The fraction of the stellar 
galaxy mass formed in-situ mi ns /m^ is highest, up to 
60%, for low mass galaxies and declines almost linearly 
(despite some scatter) with increasing galaxy mass down 
to 1=3 13% for the most massive systems in our simula- 
tions which are the central galaxies of massive groups 
(top panel of Fig. [TT]) . This trend is very simil ar to 
semi-analytical predictions (|Khochfar fc S ilk 2006a|) and 
constraints based on halo occupation models combined 
with isolated merger simulations (?). In the central panel 
of Fig. [TTJwe show the fraction of in-situ mass versus the 
concentration parameter c of the dark halo which is de- 
fined as the ratio between r2oo and r s , where r s is the 



p(r) 



$ c Per it 



{r/r s )(l+r/r s f 



(3) 



For the fit we binned the halo into 32 spherical 
shells equally spaced in /o q-in(r) between r^ n p an d 
logio(r/r2oo) = —2.5 similar to fGrossi fc Springe!! (2009) . 
We see a continuous change of the dark matter halo 
concentration. As expected from the effect of adia- 
batic contraction galaxies with significant in-situ star 
formation, i.e. more dissipation, live in more con - 
centrated halos (iBlumenthal et al.lll986t lpubinskilli"99~ 



Jesseit et al.l l2002t iDebattista et ah 
2004tlAbadi et al.ll2010HA'uger' et al.l 



2008: iGnedin et al 
20101 ). The concen 



tration of more massive halos does not increase signif- 
icantly as the matter is added pr edominantly in stellar 
form and cannot dissipate (see e.g. Uohansson et al.ll2009l 
and references therein), i.e. the adiabatic contraction ap- 
proximation cannot be applied for massive galaxies. The 
bottom panel in Fig. fTTI shows the stellar density inside 
the spherical half-mass radius versus the ratio of in-situ 
created stars of the galaxies. The two properties are cor- 
related in the sense that galaxies with a large fraction of 
accreted stars have lower central d ensities, a well know n 
property of elliptical galaxies (e.g. IBender et al.l [l992'). 

Fig. [10] gives a clue to the paradox of 'downsiz- 
ing'. The initial expectation was that in a hierarchi- 
cal universe, since more massive halos statistically are 
formed later than less massive ones, the same should 
be true of galax ies. But we know that this is not true 
observationally (jNelan et al.l I2005T) . giant ellipticals are 
older - not younger - than lower mass systems (see e.g. 
iThomas et al.l [20051 ). Our simulations give the same re- 
sult as can be seen from Fig. [151 the most massive sys- 
tems are made out of the oldest stars. The inclusion of 
galactic winds would probably lead to less efficient star 
formation at high redshifts and leave more gas for late 
in-situ star formation especially in the lower mass sys- 
tems, rendering these galaxies even younger. This would 
lead to an even steeper relation than the one shown in 
Fig. rj2J The explanation of the paradox is obvious: 
The accreted stars are typically made in smaller systems 
and these small systems are in fact made at early times 
(dashed red curves in Fig. [TO")) . Massive galaxies are 
more dominated by the accreted stars and so by z=0 they 
contain primarily old stars, although the galaxies them- 
selves are assembled late. iDe Lucia et al.l [20061 obtain 
the same result with their semi-analytic model. This way 
the expectation from hierarchical structure formation is 
satisfied. Both our simu lations and the observations of 
van DokkumlFan ((2008) agree: even at late times mas- 
sive galaxies continue to grow in mass and size. 

4. GALAXY SIZES 

The left panel of Fig. [13] shows the present-day spher- 
ical half-mass radius for the different components of our 
galaxies. The size of the in-situ component shows a very 
weak trend with galaxy mass. For the low mass galaxies 
the half-mass radii of the in-situ and the accreted stars 
are of similar size. While the in-situ component does 
not get larger than « 3kpch _1 , the half-mass radius of 
the accreted stars is strongly increasing with galaxy mass 
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Fig. 13. — Left panel: Stellar mass inside 10% of the virial radius vs. spherical half-mass radius of accreted (green diamonds), in-situ 
(black triangles) and all stars (red squares), respectively. The dashed lines show the results of a linear fit for the respective components 
( r i/2 l°g M", with a = 7.5, 6.1 and 1.1 for the accreted, total and in-situ stars, respectively). While the half-mass radius of the 
accreted stars strongly increases with mass, the half-mass radius of the in-situ formed stars shows only a weak dependence on galaxy mass. 
The mass-size relation is driven by the by the accreted stars. Right panel: This plot shows the spherical half-mass radii of the galaxies 
as a function of the ratio of accreted to in-situ created stars. The size increase of the galaxies is roughly linear dependent on this ratio 
(n/ 2 OC 0.97 * Ma.cc/Mi n3 ). 



and since the fraction of accreted stars rises with galaxy 
mass as well, the global half-mass radius of the galax- 
ies follows this trend. In our simulations the majority of 
the in-situ created stars are formed in the bulges of the 
galaxies. Stronger feedback mechanism would probably 
lead to more star formation in galactic disks resulting in 
larger radii of the in-situ component. The half- mass ra- 
dius of the accreted stars should not be affected by this. 
The right panel of Fig. [13] shows the galaxy radii versus 
the ratio of accreted to in-situ created stars. We find an 
almost linear trend. Fig. IT3l shows, that stellar accretion 
is the dominant mechanism for the size growth of mas- 
sive galaxies. For most of our systems the accretion of 
stars is significant at low redshifts, as seen in Fig. [TUJ 
especially for the high mass galaxies. Half of the total ac- 
creted stellar mass is added to the galaxies between red- 
shift one and the present day which leads to considerable 
size increase at late times. Co nsistent with this p i ctures 
are the observations from e.g. Ivan Dokkum et all (2008) 
and others that show a significant growth between red- 
shift z = 3 and z = for quiescent early type galaxies. 
A detailed analysi s of this effect will be presented sepa- 
rately. We give, in lNaab et al.l ((2009) a simple argument 
based on the virial theorem showing how late accretion 
of low mass satellites ('minor mergers') will lead to the 
rapid growth in galactic size. 

5. SUMMARY AND DISCUSSION 

We present results from 39 cosmological re-simulations 
of dark matter halos including gas and star formation 
covering a mass range of almost two orders of magnitude 
in virial mass. In the study presented here we used the 
simulations to investigate fundamental formation and as- 
sembly processes, i.e. how and when do galaxies get their 
gas and stars, and how does this influence the present day 
galaxy properties. 

We have shown that it can be useful, at a very ba- 



sic level, to distinguish between stars that are created 
inside the galaxies themselves (in-situ) and those that 
are created outside (ex-situ) and are accreted later on. 
The division into these two separate phases is quite clean 
(see Fig. [7|) with the in-situ stars typically formed closer 
than r/r v i r ~ 10 _1 , i.e. within the galaxy, and the ex- 
situ stars formed outside the galaxy at r/r V i r ~ 10 5 to 
10. Independent of galaxy mass we find that the forma- 
tion of the accreted stars peaks at redshift z ~ 4. The 
in-situ star formation as well has an early peak but ex- 
tends over a longer period of time. The ratio of stars 
that are created in-situ to the accreted stars, however 
varies strongly for galaxies of different masses. We find 
that for massive galaxies (~ 1.9 — 3.6 x 10 11 /i -1 Mq) the 
contribution of in-situ and accreted stars becomes com- 
parable early {z rj 2) and the accreted stars can account 
for up to 87% of the final stellar mass. The lower mass 
galaxies (~ 4 — 10 x l0 lo h~ 1 MQ) still can have a high 
fraction of in-situ formed stars up to 60% at the present 
day. They show a significant amount of in-situ star for- 
mation throughout the whole simulation time. The large 
difference in time when those accreted stars are actually 
formed and when they are finally assimilated by their 
host, together with the trend shown of the ratio of in- 
situ formed st ars explains the pheno menon of 'downsiz- 
ing' (see also (|De Lucia et al.l [2006) for semi-analytical 
simulations). The more massive galaxies consist mainly 
of accreted and therefore old stars leading to the depen- 
dence of mean stellar age to galaxy mass shown in Fig. 
1121 The massive galaxies in our sample assemble about 
half their mass below a redshift of z=l. This mass in- 
crease, caused by stellar accretion and merging, is not 
accompanied by significant star formation, can be a sig- 
nificant contribution to the observed increase of stellar 
mas s in the early-type galaxy popu lation since z=l (see 
e.g. iBrown et al1l2007l : iFaber et al.ll2007| ). 
We find that the accreted stars are primarily responsi- 



The Two Phases of Galaxy Formation 



11 



TABLE 2 

The assembly of stars in massive galaxies 





In-situ 


Accreted 


Epoch 


6 > 2 > 2 


3 > z > 


Baryonic mass source 


cold gas flows 


minor & major mergers 


Size of region 


r 1/2 fa 2kpc 


r 1/2 « 7kpc 


Energetics 


Dissipational 


Conservative 



ble for t he low redshift size incre ase in massive galaxies 
(see e.g. iHvde fc Bernardil 12009). When looking at the 
half-mass radii of the galaxies and the half-mass radii of 
the in-situ created and accreted components, we find that 
the half-mass radius of the in-situ created stars is only 
weakly dependent on the galaxy mass and is quite small 
(< 3kpch _1 ). This component forms at redshift z > 2 
and makes the co mpact cores of the galaxies (see e.g. 
Ivan Dokkum et ail 120081 and references therein) . The 
larger sizes of galaxies with larger mass are mainly due to 
the accreted stars creating an outer envelope with half- 
mass radii exceeding 8kpch _1 . 

Our simulations overestimate the stellar mass of the 
galaxies by roughly a factor of 2. This is probably due to 
the lack of ejcctivc and preventive feedback mechanism in 
our simulations. The stars that are accreted as well as the 
early formed in-situ stars are generated in small systems 
where winds are most effective and lead to lower star 
formation rates and therefore lower accretion rates at 
lower redshifts. The late in-situ star formation should be 
diminished by AGN feedback particularly in the massive 
systems. It will be worthwhile to investigate whether 
and how the inclusion of those processes could influence 
the presented balance of in-situ star formation to stellar 
accretion. 

The description of galaxy formation as a two phase 
process followed in a seemingly natural way from 
our detailed hydro simulations and is organized into 
a coherent scheme in Table It is not intended 
as a rival to other ways of seeing galaxy formation 



but rather as a framework within which the physical 
processes can be understood in a straightforward way. 
Early, in-situ star formation is clearly sim ilar to that 
resulting from the 'co ld flow' picture (|Dekel et al.l 
I2009al IKeres et al.l 120051) or the earlier descriptive term 
'dissipative collapse'. In fact the in-situ phase bears an 
uncanny resemblan c e to the 'monolithic collaps e ' model 
(lEggen et al.l Il962t iPartridge fc Peebles! Il967t iLarsonl 
Il969t iSearle et al.lll973t lLarsonlll975lh The late assem- 
bly phase of massive galaxies has many aspects similar 
to the 'd ry merger' paradigm inve st igated by many 
author s (IKhochfar fc Burkertl [2001 IKhochfar fe Silkl 
20061* iNaab et al.l 120061: Ivan der Wei et al.l 120091: 
Bezanson et al.l l2009t iNipoti et al.l l2009bl ) with the 
added qualification that most of the accreted stellar 
systems are low in mass compared to the final assembled 
i.e. minor mergers d ominate. There ap pears to be 
recent archaeological (jCoccato et al. I2010D and direct 
observ ational for this scenario. van Do kkum et al.l 
(2010) conclude that massive compact galaxies at z=2 
(the end of the in-situ phase) have increased their mass 
at radii r > 5kpc by a factor of « 4 since z=2 with the 
mass at smaller radii being essentially unchanged. 

We thank the referee for valuable comments on the 
paper. We are grateful to S. D. M. White, S. Faber, P. 
van Dokkum, M. Kriek , C. Conroy and M. Franx for 
helpful discussions on the manuscript. The simulation 
were partly performed at the Princeton PICSciE HPC 
center. This research was supported by the DFG cluster 
of excellence 'Origin and Structure of the Universe'. 
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